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Abstract 

This paper studies switching stabilization problems for general switched nonlinear systems. A piecewise smooth control- 
Lyapunov function (PSCLF) approach is proposed and a constructive way to design a stabilizing switching law is developed. 
The switching law is constructed via the directional derivatives of the PSCLF with a careful discussion on various technical 
issues that may occur on the nonsmooth surfaces. Sufficient conditions are derived to ensure stability of the closed-loop 
Filippov solutions including possible sliding motions. The proposed PSCLF approach contains many existing results as special 
cases and provides a unified framework to study nonlinear switching stabilization problems with a systematic consideration 
of sliding motions. Applications of the framework to switched linear systems with quadratic and piecewise quadratic control- 
Lyapunov functions are discussed and results stronger than the existing methods in the literature are obtained. Application 
to stabilization of switched nonlinear systems is illustrated through an numerical example. 
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1 Introduction 

This paper studies switching stabilization problems for 
general switched nonlinear systems in continuous time. 
The goal is to develop a constructive way to design state- 
feedback switching laws to ensure closed-loop stabilities 
including possible sliding motions. The problem is re¬ 
garded as one of the basic problems in switched sys¬ 
tems 0 and has received considerable research atten¬ 
tion in recent years El HZ], EE El, Ell ■ 

The literature on switched systems has been mainly fo¬ 
cusing on their stability analysis El El HI, EE E3 ■ De¬ 
pending on the assumptions on switching signal, these 
studies can be divided into three categories, namely, 
stability under arbitrary switching, stability under slow 
switching, and stability under state-dependent switch¬ 
ing. Various concepts and tools, such as common Lya¬ 
punov functions Hjl . dwell time |13|, and multiple Lya¬ 
punov functions [J, HU, have been proposed to study 
these problems. These analysis methods typical view 
the switching inputs as disturbances (e.g. stability un¬ 
der arbitrary or slow switching) or as being generated 
by a known switching law (e.g. stability under state- 
dependent switching). They cannot be directly applied 
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to solving stabilization problems, for which one has the 
freedom to design the switching input to stabilize the 
system. One main challenge is on the nonsmoothness 
of switching control laws, and the sliding motions that 
may occur in closed-loop trajectories. Existing stability 
analysis results provide little insight into these issues as 
they typically require the switching signal to be piece- 
wise constant [lOL fl8l| , thus eliminating the possibility of 
having sliding motions that are crucial for switching sta¬ 
bilization problems (see Examples 1 and 2 in Section 2). 

In contrast to stability analysis, switching stabiliza¬ 
tion has not been adequately studied in the literature. 
Most of the results focus on quadratic stabilization of 
switched linear systems (SLS). A switched system is 
called quadratically stabilizable if it admits a quadratic 
control-Lyapunov function (CLF) El El [25f . It has 
been shown that a SLS is quadratically stabilizable if 
there exists a stable convex combination of system ma¬ 
trices [2l M El El- Later, some works extend these 
quadratic stabilization results by employing piecewise 
quadratic CLFs. A well-known result along this direc¬ 
tion is the largest-region switching strategy proposed 
in EE El • Such a switching strategy is parameterized 
by a collection of (possibly overlapping) regions de¬ 
fined in terms of symmetric matrices and a collection of 
quadratic Lyapunov-like functions. The switching law 
synthesis problem is then formulated as an optimization 
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problem subject to bilinear matrix inequality (BMI) 
constraints. The resulting switching law only guaran¬ 
tees closed-loop stability if no sliding motion occurs. A 
different approach is proposed in a more recent study 
(0), for which the switching law is constructed using 
three classes of composite CLFs defined by taking the 
pointwise minimum, pointwise maximum, or convex 
hull of a finite number of quadratic functions. For all 
of these three classes of composite CLFs, some BMI 
conditions are derived to ensure closed-loop stability 
excluding sliding motions. It is further shown that sta¬ 
bility including sliding motions can be guaranteed for 
the case with pointwise minimum CLFs. 

The aforementioned results on switching stabilization 
have several limitations. Firstly, many of them focus on 
quadratic stabilization problems, which require the sys¬ 
tem to have a quadratic CLF. This can be conserva¬ 
tive even for SLSs as many SLSs are switching stabiliz- 
able without having a quadratic CLF [l8j]. Although sev¬ 
eral forms of piecewise quadratic CLFs 00m,m 
have been considered, the methods are often ad-hoc and 
rather specific to the particular structure of the adopted 
piecewise quadratic CLFs. Secondly, many results lack a 
systematic way to analyze sliding motions. In fact, sta¬ 
bility analysis without considering sliding motions may 
lead to incorrect conclusions about the actual stability 
behavior of the closed-loop system (see Examples 1 and 2 
in Section 2). Lastly, the majority of the existing meth¬ 
ods are developed exclusively for SLSs. Switching sta¬ 
bilization of switched nonlinear systems have not been 
adequately studied. 

In this paper, a piecewise smooth control-Lyapunov 
function (PSCLF) framework is developed to address 
the above limitations. We focus on PSCLFs for several 
reasons. Firstly, they are less restrictive than smooth 
CLFs in the sense that a switching stabilizable system 
may not admit a smooth CLF. Secondly, directional 
derivative always exists for PSCLFs, which is of cru¬ 
cial importance for constructive design of stabilizing 
switching laws. Lastly, they are general enough to cover 
a wide range of applications. In fact, most of the ex¬ 
isting nonsmooth CLFs proposed in the literature are 
special classes of PSCLFs. Examples include piecewise 
quadratic CLFs used in 1331 and the three classes of 
composite CLFs proposed in [Tdl . 

In addition to its generality, the proposed PSCLF frame¬ 
work also enables a constructive design of stabilizing 
switching laws with a systematic consideration of closed- 
loop sliding motions. In particular, we propose a gen¬ 
eral procedure to construct a switching law based on a 
given PSCLF and carefully discuss several technical is¬ 
sues when the sliding surface partly coincides with the 
nonsmooth surface of the PSCLF. We show that the re¬ 
sulting switching law always guarantees the existence of 
Filippov solutions of the closed-loop system. We also de¬ 
rive sufficient conditions to ensure closed-loop stability 


including sliding motions. Different from many existing 
works, our stability results differentiate sliding motions 
on attractive sliding surfaces from the ones on nonat- 
tractive sliding surfaces. This allows us to focus on slid¬ 
ing motions that are of practical importance, making the 
overall framework less conservative. The effectiveness of 
the framework is illustrated through several analytical 
examples for switched linear systems and a numerical 
example for a switched nonlinear system. The analyti¬ 
cal examples also lead to stronger results than the ex¬ 
isting methods on stabilization of SLSs. The proposed 
PSCLF approach, along with its stability analysis, con¬ 
tains many existing works as special cases and can be 
used as a unified framework to design stabilizing switch¬ 
ing laws with a systematic consideration of sliding mo¬ 
tions. 

It is worth mentioning that our results on switching sta¬ 
bilization can not be obtained using the classical re¬ 
sults on nonsmooth stabilization based on nonsmooth 
CLFs (e.g. P-l5l-[^ l26U29l p. This is due to the fundamen¬ 
tal differences between the two classes of problems, (i) 
Firstly, most classical nonsmooth stabilization results fo¬ 
cus on traditional nonlinear systems for which the open- 
loop vector field f(x,u) is continuous in ( x,u ). These 
results are not directly applicable to switched systems 
whose vector field is not continuous with respect to the 
switching control, (ii) Secondly, while we focus on Filip¬ 
pov solutions due to its crucial importance for switched 
systems, many works on nonsmooth stabilization adopt 
sample-and-hold solutions (e.g. d-Q) or Caretheordory 
solutions (e.g. 0 )- The difference in the adopted solu¬ 
tion notions will significantly change the nature of the 
problems, making previous methods not directly appli¬ 
cable to our problem, (iii) Lastly, many classical results 
merely provide sufficient and/or necessary conditions for 
nonsmooth stabilizability. For example, the classical re¬ 
sults on nonsmooth CLF (e.g. I^-Q) show that asymp¬ 
totic controllability of a nonlinear system is equivalent 
to the stabilizability by a measurable feedback law in 
the sample-and-hold sense. However, there is little indi¬ 
cation on how to find such a stabilizing law in a con¬ 
structive and computationally tractable way. In con¬ 
trast, this paper offers a simple and constructive way to 
design the switching law based on a given PSCLF and its 
corresponding partitions. Due to the above differences, 
switching stabilization cannot be viewed as a special case 
of the classical nonsmooth stabilization problems, and 
the proposed PSCLF framework represents a novel con¬ 
tribution. 

The rest of the paper is organized as follows: In Section 
2, we use two examples to illustrate the importance of 
sliding motions and motivate the problem formulation. 
In Section 3, we introduce the PSCLF framework and de¬ 
rive two PSCLF theorems to ensure closed-loop stability 
excluding and including sliding motions, respectively. In 
addition, we discuss in detail two important special cases 
with smooth and pointwise minimum CLFs. In Section 
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4, we use the proposed framework to study stabilization 
of SLSs under three types of CLFs and provide a nu¬ 
merical example for stabilization of switched nonlinear 
systems. 

Notations: Let R + be the set of nonnegative real num¬ 
bers, R” be the n-dimensional Euclidean space. Denoted 
by dM and M the boundary and the closure of a set 
M C R n , respectively. For any positive integer m , de¬ 
noted by N m = {1, ■ ■ • , m} the set of positive integers 
that are less than or equal to m. Denoted by | -1 the car¬ 
dinality of a given set, and || • || the Euclidean norm of 
a given vector or matrix. Let p be Lebesgue measure. 
For any x G R” and e > 0, define J\f(x\ e) = {z G R n : 
lk-^ll < r}- 

2 Problem Formulation 

In this paper, we consider the following continuous-time 
switched nonlinear system 

x(t) = f a ( t) (x(t)), ff(i)eS = {l,-,M}, (1) 

where x(t) G R" denotes the continuous state of the sys¬ 
tem, fi : R” —>• R” denotes the vector field of subsystem 
iG Q, and aft) denotes the switching control signal that 
determines the active subsystem at time t G R+. We as¬ 
sume that the vector field of each subsystem fj , i G Q is 
locally Lipschitz continuous on R™. In addition, the ori¬ 
gin is assumed to be a common equilibrium for all the 
subsystems {fi} l( zQ- 

Assume that the state x(t) is available at all time t G R+, 
and the switching control is determined through a state- 
feedback switching law v : R” —> Q. The corresponding 
closed-loop system can be written as 

= fv{x{t))(. x (t)) ■ (2) 

For simplicity, the closed-loop vector field under 
a switching law v will be denoted by i.e., 

f\ v ](x) = fv( x ){x), \/x G R”. Although each subsystem 
vector field fj is continuous, a switching law v may 
introduce discontinuities in the closed-loop vector field 
/[„]. In general, the differential equation in (2) may not 
have a classical or Caratheodory solution 0. In this pa¬ 
per, we adopt the Filippov solution notion [TTJ to study 
the switching stabilization problem. 

Definition 1 (Filippov Set-Valued Map f^j) For any 
vector field X : R" —> R", the corresponding Filippov 
set-valued map F[X] : R n —> *8(R") is defined as 

F[X](x)=p| p| cd{X(J\f(x\ <5)\5')} , x G R", 
<5>0 A j(S)=0 


where f8(R”) denotes the collection of subsets of R™, co 
denotes convex closure and p denotes Lebesgue measure. 

Definition 2 (Filippov Solution J^j) A Filippov solution 
to a differential equation x(t) = X(x(t)) over [0,ti] with 
ti > 0 is an absolutely continuous map x : [0,1 1 ] —> R ra 
that satisfies the differential inclusion x{t) G F[X](x(t)) 
for almost all t G [0,ii]. 

According to Definition 1, we can exclude an arbitrary 
set of measure zero around x when computing F[X](a;). 
This implies that two vector fields that differ on a set of 
measure zero will lead to the same Filippov set-valued 
map, and hence the same Filippov solution. In addi¬ 
tion, it can be verified that if the vector field A' is con¬ 
tinuous, then the Filippov solution to x(t) = X(x(t)) 
coincides with the classical solution. If A' is discontin¬ 
uous, then the Filippov solution exists as long as the 
map X : R" —► R n is measurable and locally essentially 
bounded EH- Since we assume that the vector field of 
each subsystem is continuous, it can be easily verified 
that the Filippov solution to the closed-loop system (2) 
exists whenever the switching law v : R” —> Q is mea¬ 
surable. However, switching laws may take unappealing 
forms if we merely require them to be measurable. For 
the purpose of systematic analysis on sliding motions, 
we focus on a class of switching laws defined as follows. 

Definition 3 (Admissible Switching Law) A switching 
law v : R" —> Q is called admissible if there exists a col¬ 
lection of disjoint and open sets {Dj}i g g C R" satisfying 
U ieQVi = R” such that v(x) = i, Vx G Vj. 

The sets {Vi} i& Q and {dVi} i( zQ in Definition 3 will be 
referred to as the switching partitions and the switch¬ 
ing boundaries, respectively. Let dT> = Uj g g<92};. In 
the above definition, each switching partition T>i is not 
required to be connected. In addition, since {Tk}j g g 
are full-dimensional open subsets of R", the switching 
boundaries {92>i}j g g are of measure zero. As a result, 
an admissible switching law v is piecewise constant. It 
is easy to see that an admissible switching law ensures 
the closed-loop vector field /[„] to be measurable and lo¬ 
cally essentially bounded. Thus, a Filippov solution to 
system (2) always exists. 

We denote x(-\ z, v) : R + —>• R n as a Filippov solution to 
the closed-loop system (2) under an admissible switching 
law v with initial state z G R”. Sliding motion may 
occur along a Filippov solution [32j. More precisely, if 
there exists a nontrivial time interval (^ 1 ,^ 2 ) such that 
x(t; z, v) f[w](x(t; z , v)) for almost all t G (ti,^ 2 ), then 
this part of the state trajectory is called a sliding motion 
where the velocity of the state z, v) can take any 
value within the convex set F[fy v (\(x(t\ z, v)). 

Since the closed-loop vector field /[„] is continuous inside 
each of the switching partition T>i, a sliding motion of 
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/[„] can only occur on the switching boundaries dV. For 
each x £ dV, let X sm (x) be the set of indices of all the 
subsystems that may possibly be involved in a sliding 
motion, i.e., 


X sm (x) ±{i£Q:xe dVi}. (3) 

During sliding motion, the velocity satisfies x = 
Sigi m (x) a i{ x )fi( x ) for some convex combination co¬ 
efficients {«z(a;)} ie x sm (x) with Ei e x am ( x ) a *( x ) = 1 - 
These coefficients can be determined if the sliding 
surface equation is known |To|, [3l|. In addition, we 
define If m (x) as the set of indices of all the ac¬ 
tive subsystems involved in the sliding motion, i.e., 
X ’%m(x ) = {i £ x sm {x) : cti(x) ± 0 where a; = 
ism{x) a i( x ) = !}• Obviously, 

lsm( x ) C X sm (x), Vx £ dV. 

Existing works in the literature often discuss stability 
properties by excluding sliding motions. However, this 
may lead to an incorrect conclusion about the actual 
stability behavior as illustrated in the following example. 

Example 1 (Unstable Sliding Motion I) Consider the 
following piecewise linear system: 



—3 -l' 


12 2 


" 


-3 1 


-12 2 


x (t), 


x(t), 


if x\(t) > 0 
if x\ (t) < 0 


(4) 


The closed-loop vector field is discontinuous on the sur¬ 
face {x £ R 2 : x\ = 0}, resulting in a sliding motion. 
It is easy to see that both subsystems are stable and the 
closed-loop system will also be stable if no sliding mo¬ 
tion occurs. However, the phase portrait in Fig. 1(a) in¬ 
dicates that the upper part of the switchin g s urface (i.e. 
{x £ R 2 : xi = 0, X2 > 0}; is attractivAAand all the 
Filippov solutions will involve sliding motions regardless 
of the initial state locations. Fig. l(b) shows two trajec¬ 
tories starting from xq = [0.5,—2]” simulated with two 
hysteresis bands of different sizes. Both of them involve 
an unstable sliding motion that grows unbounded along 
the X 2 coordinate. 

Example 1 indicates that excluding sliding motions may 
lead to an incorrect stability conclusion for switched sys¬ 
tems. Therefore, sliding motions should be carefully con¬ 
sidered during the switching law design. On the other 
hand, simply requiring stability for all sliding motions 


1 Roughly speaking, a sliding surface is attractive if all the 

trajectories starting from a sufficiently close neighborhood 
will converge to the surface. When the surface is differen- 
tiable, simple conditions are available to check its attractive¬ 
ness [321], 


can be overly restrictive for many applications as illus¬ 
trated in the next example. 

Example 2 (Unstable Sliding Motion II) Consider the 
following piecewise linear system: 



’-3 -1 


12 2 


" 


-3 1 


-12 2 


x(t), 


x(t), 


if xi(t) < 0, X 2 (t) > 0 
or Xi(t) > 0, X 2 (t) < 0 

if Xi(t) > 0, X2 (t) > 0 
or xi (t) < 0, X 2 (t) < 0 


(5) 


As illustrated in Fig. 1(c), the closed-loop vector field is 
discontinuous on the surface Si U S 2 where Si = {x £ 
R 2 : X\ = 0} and S2 = {x £ M 2 : X2 = 0}. Accord¬ 
ing to the phase portrait in Fig. 1(c), if the initial state 
Xq £ R 2 \5i, the trajectory will eventually reach S 2 fol¬ 
lowed by a stable sliding motion converging to the origin 
along S2. If the initial state Xq £ Si, the Filippov solu¬ 
tions are nonunique. Theoretically speaking, both stable 
and unstable Filippov solutions may exist (see Fig. 1(c) 
for illustration). However, the sliding surface Si, corre¬ 
sponding to the unstable sliding motion, is not attrac¬ 
tive. As a result, any numerical simulation with nontriv¬ 
ial hysteresis bands can only produce stable trajectories 
as shown in Fig. 1(d). In other words, the ideal unsta¬ 
ble sliding motion (green dashed trajectory in Fig. 1(c)) 
cannot occur in real systems. 

In Example 2, all the Filippov solutions converge to the 
origin except the sliding motions on Si. Since the slid¬ 
ing surface Si is not attractive, sliding motions on this 
surface cannot sustain in reality and should not be con¬ 
sidered in the analysis. Motivated by the above two ex¬ 
amples, we will define two sets of Filippov solutions with 
one excluding sliding motions and the other including 
the sliding motions on attractive sliding surfaces (Exam¬ 
ple 1) but not the ones on nonattractive sliding surfaces 
(Example 2). 

Definition 4 Given an admissible switching law v and 
a positive constant e > 0, the set of Filippov solutions 
excluding sliding motions is defined as 

X 0 {e,v) = {x(-;z,u) : p|| < e, 

x(t\z,v) = f [u] (x(t-,z,iy)),fort £ M+ a.e.} 


Before introducing the set of Filippov solutions includ¬ 
ing sliding motions, we first establish the relation be¬ 
tween the switching partitions and the vector fields of 
the active subsystems involved in the sliding motion by 
the following regularity condition. 

Definition 5 Assume x £ dT>, x is said to be regular if 
Vi £ Xf m (x), there exists ei > 0 such that x—6fi(x) £ Vi, 
V<5 £ (0, ej). 
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(a) 




(b) 



(c) (d) 

Fig. 1. (a): The phase portrait of system (4), (b): The closed-loop trajectories of system (4) starting from xo = [0.5, — 2] T , (c): 
The closed-loop trajectories of system (5) starting from xo = [0, —0.5] T , (d): The closed-loop trajectories of system (5) with 
hysteresis band A = 0.1. 


Roughly speaking, if x £ dT>i is regular, then for any 
subsystem vector field fi ( x ) that contributes nontrivially 
to the velocity x (i.e. a,;(x) ^ 0), we can conclude that 
fi{x) must point outwards T>i. Note that if the switch¬ 
ing boundaries dT> are continuously differentiable, reg¬ 
ularity condition can also be stated as V* £ Tf m (x), 
d(j)i(x) T ■ fi{x) > 0 where (j>i{x) is a differentiable func¬ 
tion such that Dj = (i t I" : (fi{x) < 0}. Now we are 
ready to introduce the set of Filippov solutions that con¬ 
tains not only all the solutions in Xg but also the type 
of sliding motions in Example 1. 

Definition 6 Given an admissible switching law v and 
a positive constant e > 0, the set of Filippov solutions 
including sliding motions is defined as 

X sm (e,u) = {x(-,z, v) : j|^]| < e,if3(t 1 ,t 2 ),ti < t 2 
s.t. x(t-z,n) ^ f{ v }{x{t\z,v)) fort £ (ti,t 2 ) a.e., 
then x(t\ z , v) is regular \/t £ {ti,t 2 )} 


We proceed to define asymptotic stability of the closed- 
loop system excluding and including sliding motions, re¬ 
spectively. 


Definition 7 System. (2) under an admissible switching 
law v is asymptotically stable excluding sliding motions 
if (i) for any e > 0, there exists S > 0 such that for 
any Filippov solution x (•; z, v) £ Xg (<5, v), ||x(f; z, v)\\ < 
e, Vi £ R+. (ii) lim^oo x(t; z, v) = 0. 

Definition 8 System (2) under an admissible switching 
law v is asymptotically stable including sliding motions 
if (i) for any e > 0, there exists S > 0 such that for any 
Filippov solution x{-\z,v) £ X sm (5,v), ||x(t; z, i^)|| < 
e, Vt £ R+. (ii) lim^oo x(t; z, v) = 0. 

With these definitions, we are ready to formally state 
the switching stabilization problem to be studied in this 
paper. 

Problem 1 (Switching Stabilization Problem) Find an 
admissible switching law v under which the closed-loop 
system (2) is asymptotically stable including sliding mo¬ 
tions. 

While there are many ways to study the switching stabi¬ 
lization problem, we are particularly interested in meth¬ 
ods that (i) can lead to constructive ways to find a sta¬ 
bilizing switching law, (ii) allow for systematic stabil- 
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ity analysis for sliding motions, and at the same time 
(iii) are less conservative and more general than existing 
methods. In the rest of this paper, we will first develop a 
control-Lyapunov function framework to accomplish the 
aforementioned objectives and then discuss its connec¬ 
tions to the literature and its applications to switched 
nonlinear systems. 

3 A Piecewise Smooth Control-Lyapunov Func¬ 
tion Framework 

The goal of this section is to develop a unified control- 
Lyapunov function (CLF) framework to solve the sta¬ 
bilization problems for switched nonlinear systems. We 
focus on piecewise smooth control-Lyapunov functions 
(PSCLFs), which allows for constructive design of stabi¬ 
lizing switching laws. To present the framework, we will 
first introduce the concept of PSCLFs, and then develop 
a strategy to construct switching laws using PSCLFs. 
Sufficient conditions will be derived to ensure closed-loop 
asymptotic stability including sliding motions. We will 
also discuss two important special cases of our PSCLF 
framework, namely, smooth CLFs and pointwise mini¬ 
mum CLFs, for which we will show that stable sliding 
motions are always guaranteed without any additional 
requirement on the sliding surfaces. 

3.1 Piecewise Smooth Control-Lyapunov Function 

Nonsmooth CLFs have been well studied for stabiliza¬ 
tion of traditional nonlinear systems [••y affil [2S| . Their re¬ 
lation to our switching stabilization problems has been 
discussed in the introduction section. This subsection 
will introduce a particular class of nonsmooth CLFs, 
namely, piecewise smooth CLFs, for switching stabiliza¬ 
tion problems. We first give a formal definition for piece- 
wise smooth functions. 

Definition 9 (Piecewise Smooth Function) A function 
g : R n — > R is called piecewise smooth if it is continuous 
and there exists a finite collection of disjoint and open 
sets, ffi,..., Cl m C R" such that (i) Uj e pj m f2j = R n , (ii) 
g is continuously differentiable onQj, Vj £ N m , and (iii) 
dflj is a differential manifold for each j £ N m . 

Remark 1 Piecewise smooth functions can be defined in 
different ways IBM]- In our definition, each partition 
f lj is not required to be connected. In addition, since each 
Qj is a full-dimensional open set in 1", its boundary dflj 
must have measure zero, i.e. nfdClj) = 0. 

An important property of piecewise smooth function is 
the existence of directional derivative anywhere in the 
state space as stated in the following lemma. 


tion 9 is directionally differentiable on M™, i.e., the limit 
Dg{ cc; 77 ) = lim^o j(g{x + Srj) — g(x)) exists,Vx,q £ R". 

Control-Lyapunov function is a useful tool to study sta¬ 
bilization problems. To be less conservative, we will fo¬ 
cus on piecewise smooth control-Lyapunov functions for 
general switched nonlinear systems ( 1 ) defined as fol¬ 
lows. 

Definition 10 (Piecewise Smooth Control-Lyapunov 
Function (PSCLF)) A piecewise smooth function 
V : R n —> R + is called a PSCLF if there exists an¬ 
other continuous function W : R™ —> R + such that the 
following conditions hold: 

V(x) > 0, W(x) > 0, Vx yf 0, V (0) = 0 ( 6 ) 

C /3 = {x : V(x) < /?} is bounded for each f3 (7) 
min DV(x; fi(x)) <—W(x), Vx £ R" ( 8 ) 

In addition, the pair of functions (V, W ) is called a CLF 
pair. 

Remark 2 The PSCLF defined above can be viewed as 
a special class of nonsmooth CLFs fH /. Roughly speak¬ 
ing, the existence of such a function implies the stabi- 
lizability of the system in the sample-and-hold sense W- 
The main complication for switching stabilization lies in 
how to constructively find a stabilizing switching law in 
the Filippov sense and how to guarantee the closed-loop 
stability including the sliding motions on the attractive 
sliding surfaces. 

We refer to ( 6 ) as the positive definite condition, refer 
to (7) as the radially unbounded condition, and refer to 
( 8 ) as the decreasing condition. Note that the decreas¬ 
ing condition is given in terms of directional derivative, 
which is well-defined for piecewise smooth functions ac¬ 
cording to Lemma 1. In the rest of this section, we will 
first develop a constructive way to design a switching law 
using a given PSCLF, and then carefully study various 
stability properties of the closed-loop system. 

3.2 Switching Law Construction 

If system (1) admits a PSCLF V, then it can be used 
to construct a stabilizing switching law. The main idea 
is to select the subsystem along whose vector field the 
PSCLF decreases at the fastest rate. Naturally, one may 
want to directly construct the switching law as 

2 >(x) = argmin DV(x; fi(x)), Vx £ R". (9) 

ieQ 

The switching law u : R” —>■ 2® defined above is set¬ 
valued and the part of the state space corresponding to 
multiple minimizers (i.e. X s = {x £ R" : |z>(x)| > 1}) 


Lemma 1 (Directional Derivative \3(\ l) A piecewise 
smooth function g : R™ —> R as defined in Defini¬ 
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is where sliding motions may possibly occur. To under¬ 
stand the stability behaviors of all the closed-loop Filip¬ 
pov solutions, there are several technical issues requiring 
special attention. First, the set X s may not have mea¬ 
sure zero, which complicates the derivation and analysis 
of the Filippov solutions of the closed-loop system. In 
addition, even if X s is of measure zero, it may intersect 
the nonsmooth surface of V, introducing additional chal¬ 
lenges in analyzing the sliding motions. Some aspects of 
these technical issues regarding the switching law D are 
further discussed in Example 3 in Section 3.3. To better 
address these issues, we will develop a slightly different 
method to construct the switching law, which can facil¬ 
itate our discussion on Filippov solutions and their sta¬ 
bilities. 

We first introduce some notations. Suppose we are given 
a general PSCLF V : R n —> R + with partition sets 
{Oj jJLi ■ According to Definition 9, each partition Llj is 
an open set in R™ with boundary dfl 3 and closure 0 ? -, 
j £ N m . The union of all the partition boundaries will 
be denoted by = UjeN m dil 3 . Let J(x) be the set of 
partitions whose closure x belongs to, i.e., 

J(x ) = {j £ N m : x £ Dj}, x £ R”. (10) 

J{x) is a singleton if x lies inside one of the_partitions, 
and is set-valued if x £ <9f 1. Denote by V 3 : Ll 3 —>• R the 
restriction of the function V to the closure f\ 3 C R”. 
By Definition 9, Vj is continuously differentiable on f l 3 . 
Therefore, we have 

DVj(x\r /) =< XVj(x),r] >,Vx £ % and p £ R", (11) 

where VV(x) denotes the gradient of V at x. For 
any boundary point x £ dCt 3 , we define Wj(x) = 
Ym\ xe o,. x ^%S7Vj(x). Since V 3 is continuously differen¬ 
tiable on f lj , we have 


(3) Construct the switching law as: 


v(x) 


i, if 3 i £ Q , s.t. x £ Vi 
min{* £ Q : x £ dV ,;}, if x £ dV 


(14) 


The above procedure will be represented by an operator, 
denoted by S , which maps a PSCLF V to a particular 
switching law v = <S[V]. Due to the continuities of all 
the vector fields and VV 3 , the set V 3l defined in (13) is 
an open set in R" for each j £ N m and i £ Q. Hence, 
the set Vi is also open for each i £ Q. In addition, the 
boundaries dVi, i £ Q, and their union dV are all of 
measure zero. In summary, the sets {Thj-igg are open 
and disjoint. Therefore, v defined in (14) is an admissi¬ 
ble switching law with switching partitions {Th}j g g and 
switching boundaries {(9X\}i e g. 

It is also easy to see that the switching law v defined 
in (14) assigns each state x £ R ra to a unique subsys¬ 
tem. It is worth mentioning that there are many differ¬ 
ent ways to design the switching laws on the boundary 
set dV. All of them will lead to the same closed-loop 
Filippov set-valued map as dV is of measure zero. The 
particular construction given in (14) ensures that v{x) 
at any boundary point x £ dV is a continuous extension 
of v on one of the partitions whose boundaries touch x. 

Next, we evaluate the closed-loop system (2) under the 
switching law v defined in (14). Clearly, the closed-loop 
vector field /[„] is measurable and locally essentially 
bounded. Thus, Filippov solutions always exist. In ad¬ 
dition, /[„] is continuous inside each switching partition 
Vi and piecewise continuous around any boundary point 
x £ dV, thus the Filippov set-valued map is given by @: 


F [f[u]\ (z) 


{fi{x)}, if 3* £ Q, s.t. x £ Vi 
W{fi{x) :i £ Q,x £ Vi) , o/w 


Vj(z) - Vj(x) =< VV 3 {x),z-x> +o{\\z - x||), 

— o (ii z — x in 

Mx,z £ Qj,j £ N m , where lim —r— --— = 0. (12) 

Hz — cc|| 


The rest of this section will analyze and discuss various 
stability properties of the closed-loop Filippov solutions 
with Filippov set-valued map F [/[„]] defined in (15). 


With these notations, the switching law can be con¬ 
structed as follows: 

(1) For each j £ N m , i £ Q, define 

V jii ±{x£n j :(VV j (x),f i ( x))< 

min {VVj(x)Jq(x))\] (13) 

qeQ,q^i ) 

(2) For each i £ Q, define V, = Vj ti and dV = 

Ui^QdVf, 


3.3 Piecewise Smooth Control-Lyapunov Function 
Theorems 

In this subsection, we will establish conditions to en¬ 
sure stability for the closed-loop system. The main re¬ 
sult consists of two parts. First, we will show that if V 
is a PSCLF, then the closed-loop system (2) under the 
switching law v = <S[V] is asymptotically stable pro¬ 
vided there is no sliding motion. Second, we will derive 
an additional condition for the PSCLF V to guarantee 
closed-loop stability including sliding motions as defined 
in Definition 8. 
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Throughout this subsection, we will assume (V, W) is a 
CLF pair, where W : R™ —> R + is a continuous non¬ 
negative function, and V : R" —»• R.+ is a PSCLF with 
partitions and nonsmooth boundaries <9f 1 = 

U In addition, the switching law v = <S[V] is 
generated from V as described in the previous subsec¬ 
tion with switching partitions {Pi}ieQ and switching 
boundaries dV = Ui 6 gc®i. 

Roughly speaking, the key of our stability analysis is to 
ensure the PSCLF V decreases along any closed-loop 
trajectory x(t). If there is no sliding motion, we have 
x(t) = fu(x(t)){x{t)) , for almost all t £ R + . In this case, 
we need to check the directional derivative along the 
closed-loop vector field, i.e. DV(x; f v ( x )(x)), x £ R". 
On the other hand, if there is a sliding motion, then 
there exists a nontrivial time interval {t\,t 2 ) such that 
x{t) = J2iex em ( x (t)) ai(x (£))/, (x(t)), for almost all t £ 
(t\,t 2 ). In this case, merely looking at DV(x, f v ( x )( x )) is 
no longer enough; we need to guarantee that V decreases 
along the switching boundaries dD. In other words, we 
need to further check DV (x\^2, i( - x ( x ) a i( x )f‘( x )) f° r 
x £ dT>. Before presenting the main stability results, we 
first derive some useful properties for these two types of 
directional derivatives. 

Lemma 2 For any point not on the nonsmooth bound¬ 
aries, i.e. x £ dfl, we have 

D V{x]f„( x ){x)) = min DV(x;fi(x)) < -W{x). 


PROOF. Assume x £ Llj for some j £ N m . Accord¬ 
ing to the property of V in (11) and the construction 
of v in (13), DV{x\f v ( x) {x)) = (VV ){x), f v(<x) {x)) < 
(VVj(x),fi(x)) = DVj(x; fi(x)), Vi £ Q. In other 
words, DV(x-, f v ( x )(x)) = min ieC DVj(x; fi(x)). It fol¬ 
lows from the decreasing condition of V in (8) that 
min i e QDVj(x; fi(x)) < —W(x), Vx £ f lj. Therefore, 
we conclude that DV(x; f„( x )(x)) < —W(x), Vx £ C2j 
for some j £ N m . □ 

As mentioned at the beginning of Section 3.2, one of the 
reasons to use the switching law v defined in (14) in¬ 
stead of the more natural choice v defined in (9) is that 
the latter one introduces additional challenges in ana¬ 
lyzing the sliding motions on the nonsmooth surface. 
In general, the two switching laws give different results 
on the nonsmooth surface. In particular, on the nons¬ 
mooth surface, the directional derivative along closed- 
loop trajectory under v differs from that under v, i.e., 
DV(x\f v ( x )(x)) DV(x\ fo( x )(x)), and the set of sub¬ 
systems that achieve the minimal decreasing rate under 
v differs from the one under v , i.e., I sm (x) ^ v{x). To 
better illustrate the subtle differences between the two 
switching laws on the nonsmooth surface, we give the 
following example. 



Fig. 2. Subsystem vector fields on the nonsmooth surface of 
Example 3 

Example 3 Consider a switched linear system, x = AiX, 
i £ Q = {1, ••• ,4} and given a PSCLF V = 
x T PjX,\/x £ f lj, j £ N2, where 
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The nonsmooth surface of V is given by {x £ R 2 : 
x T {P 1 —P 2 )x = 0} = {x £ R 2 : Xi — X 2 = 0}U{x £ R 2 : 
Xi+3x2 = 0}, which consists of two lines that intersect at 
the origin. In this example, DVj(x ; fi(x)), j £ N 2 , i £ Q 
takes the form of DVj(x\ fi{x)) = x T (AjPj + PjAi)x. 
We pick a point of the form z = [c, c) T £ dfl and ana¬ 
lyze its directional derivatives. Partitions {fb,-}j g N 2 and 
subsystem vector fields at z are shown in Fig. 2. The di¬ 
rectional derivatives at z can be easily calculated as 

on 

DV\{z\ fi(z)) = —20c 2 , DVi(z; / 2 (z)) = - y c 2 , 

DVx{z-,h{z)) = -16c 2 , Wi(z;/ 4 (z)) = 4c 2 , ( 16 ) 

DV 2 (z; fi(z)) = -12 c 2 ,DV 2 {z;f 2 {z)) = -20c 2 , 
DV 2 {z-h{z )) = —12c 2 , DV 2 {z] fi{z)) = 0. 

The above results indicate the following relations among 
the directional derivatives: 

DV 1 (z-J 1 (z)) = DV 2 (z-J 2 (z)) 

DVi(z; fi(z)) < DV\{z\ fi{z)), Vi ^ 1 (17) 

DV 2 (z-, f 2 (z)) < DV 2 {z; fi{z)), Vi ^ 2 

According to the switching law v defined in (If), the min¬ 
imal directional derivative at z is 

min min DVj(z-, fi(z)) = —20c 2 , 

j€N 2 i£Q 

which is actually the minimum of the eight terms in (16), 
and the corresponding set of subsystems that achieve this 
minimum is I sm (z) = {1,2} (solid and dashed blue ar¬ 
rows in Fig. 2). However, if using the commonly used 
















switching law v defined in (9), the minimal directional 
derivative at z is 

min DV(z;fi(z)) = mm{DVi{z; f 2 (z)), 
DV 1 (z-,f 4 (z)),DV 2 (z-J 1 (z)),DV 2 (z-,f 3 (z))} = -12c 2 

where the first equality follows from the definition of di¬ 
rectional derivative. The set of subsystems that achieve 
this minimum is v(z) = {1,3} (dotted and dashed or¬ 
ange arrows in Fig. 2). In summary, the switching law 
v defined in (If) gives a different result from the switch¬ 
ing law v on the nonsmooth surface where PSCLF V de¬ 
creases faster under v than under v. 

In general, DV{x; f u ( x )(x)) ± mini DV(x; fi(x)) on 
dfl as illustrated in Example 3. The expression for 
DV{x\ f v (x) ( x )) on the nonsmooth boundary dll can 
be quite involved. To address this issue, we introduce a 
set-valued map defined below. 

Definition 11 For any j £ N m , the set-valued map 
Mj : dflj —» 2® is defined as 

Mj(x) = {i£ Q: 3{^ fe } fe >o £ Clj n V, 

such that lim z k = x},x G dClj. (18) 

k—> oo 



Fig. 3. Illustrating example for Mj(x) 

Note that in the above definition Vi and flj are both 
open sets, so is their intersection. The above definition 
indicates that x £ dflj is a limit point of Llj 0 Vi for all 
i £ Mj (x). Roughly speaking, Mj (x) contains the set of 
subsystem indices that are chosen by the switching law 
v within an arbitrarily small neighborhood Al[x\ e) fl Clj , 
for all sufficiently small e > 0. For example, consider the 
partitions shown in Fig. 3, where Cli is the same as V\ 
and Cl 2 is the union of V 2 and V 3 . For this particular 
example, we have Mi(xi) = {1}, M 2 {x\) = {2,3}, 
Mi(x 2 ) = {1}, M 2 (x 2 ) = {3}. 

The following properties of Mj are helpful for deriving 
the main stability results. 

Lemma 3 Let x £ dCl and let J{x) be defined in (10). 


i) If x £ Vi for some i £ Q, then Mj(x) = {i}, 
Vj G J(x). 

ii) If x G dV, then\Jj e jr x 'jMj(x) =T sm (x). 

PROOF. Part i) of this lemma follows immediately 
from Definition 11 and the fact that the sets {Vi}i G Q 
are open and mutually disjoint. For part ii), it suf¬ 
fices to show that if i £ X Bm (x), then i G Mj(x) 
for some j G J(x). We thus fix an arbitrary sub¬ 
system index i G X sm (x). Clearly, x G dVi and 
3<5i > 0, such that Af(x, e) fl V t is a nonempty open 
set Ve £ (0,(ii). In addition, according to the defini¬ 
tion of J{x), we know there must exist S 2 > 0 such 
that AT(x,e) = Uje j( x ) (flj fl Af{x, e)), Ve £ (0, S 2 ). 
Therefore, 3j G J(x) such that Af(x, e) fl V t fl Cl-- 0, 
which in turn implies AI(x, e) (~l Vi fl fl-- is a nonempty 
open set for all e £ (0,do), where 6 q = min{i5i, <52}- Let 
e k = min{-^, d 0 } and pick any z k G Af(x, e k ) D V t n H-. 
Clearly, z k —> x, as k —> 00 . Therefore, i G Mj(x), 
which completes the proof. □ 


Some properties of directional derivatives on the nons¬ 
mooth surface can be revealed with the help of the set¬ 
valued map Mj defined in Definition 11. 

Lemma 4 For any j G N m , x G dClj, and q G Mj{x), 
we have 

i) DVj(x] f q {x)) = mm ieQ DVj(x; fi(x)); 

ii) DVj{x\ f q (x)) < -W(x). 


PROOF. We fix a; £ dflj for some j G N m and 
q G Mj(x). By Definition 11, 3{^fc}fe> 0 £ flj D V q such 
that linifc-^oo z k = x. Based on the definition of v in 
(14), u(zk) = q, which implies DVj(z k \ f v ( Zk ){zk)) = 
DVj(zk\ f q (zk))- Apply Lemma 2 to {z k }k >0 G Clj, 
DVj(z k ; fu{z k ){z k )) = mm i€ QDVj(zk-, fi{zk)) < 
—W{z k )• From above, we have the following relation. 

DVj(z k ; f q (z k )) = min DVj(z k ; fi{z k )) < ~W(z k ) (19) 
ieQ 

Due to the continuity of DVj(-, f q (-)) for all j £ N m , 
q G Q and the continuity of min, gi for finite number 
of continuous function gi, the equality part in (19) still 
holds in the limit as k —> 00 , which completes the proof 
of part i). Part ii) directly follows from the inequality 
part in (19) and the continuity of W(-). □ 

Although it is not generally true for x G dll that 
DV(x;f„( x )(x)) = min ieQ DV (x; fi(x)) (see Example 
3), a similar conclusion can still be obtained if we look 
at the part of nonsmooth surface that does not contain 
switching boundaries. 
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Lemma 5 For any x £ dfl\dT>, we have 
DV{x,f v ( x ){x)) < —W(x). 


PROOF. Assume x £ V q for some q £ Q. According 
to part i) of Lemma 3, q £ Mj(x), Vj £ J(x). By the 
definition of J(x), x £ dflj, Vj £ J{x). By part ii) of 
Lemma 4, DVj(x; f q {x)) < —W(x), Vj £ J(x). By the 
definition of directional derivative, for any g £ R" and 
x £ dfl, DV(x\g) = DVj*(x\g ) for some j* £ J(x). 
The desired result then follows from DV( x; f v (x){x)) = 
DVj* (x; f v ( x )(x)) for some j* £ J(x). □ 


With the above lemmas, we are now ready to state our 
first stability result. 

Theorem 1 (PSCLF Theorem I) The closed-loop sys¬ 
tem (2) under the switching law v = <S[V] is globally 
asymptotically stable excluding sliding motions. 


PROOF. Let x(t) be the closed-loop trajectory start¬ 
ing from an arbitrary initial state Xq £ R" under the 
switching law v = <S[V]. Since there is no sliding motion, 
we have x(t) dT> for almost all t £ R+. Thus, we have 
j-t v {x{t)) = DV (x{t);f Hx{t ))(x(t))) < -W{x(t)), for 
almost all t £ R+, where the last inequality is due to 
Lemma 2 and Lemma 5. The rest of the proof follows di¬ 
rectly from the classical Lyapunov theorem proof ( fl6l . 
Theorem 4.1]) by replacing the Lie derivative with the 
directional derivative DV (x; x). □ 


The next goal is to analyze the stability of sliding mo¬ 
tions that may occur along a closed-loop trajectory. We 
first introduce another technical lemma regarding the 
directional derivative along a tangent direction of the 
nonsmooth boundary. 

Lemma 6 For any x £ dfl and j, j' £ N m satisfy¬ 
ing T x (dflj) n T x (dClj’) / 0, we have DV{x\rj) = 
DVj (x; rj) = DVj<{x\rj), Vg £ T x {dDj) n T x {dDj>). 


PROOF. We pick an arbitrary point x £ dtt and let 
j,j' £ N m satisfy T x (dLlj) (~l Tx{d£lj<) ^ 0. It is easy 
to see that j,f £ J(x). We then pick an arbitrary vec¬ 
tor g £ T x {d£lj) n%{dn r ). According to the definition 
of tangent vector, for any S > 0, there exists a func¬ 
tion h : R —>• R™_with lim^o ||/i(<5)|| = 0 such that 
x + (g + h(S))S £ Clj. By the definition of J(x) in (10), 
we have x £ f lj. It follows from (12) that Vj(x + (g + 
h(S))S)-Vj(x) = (VV J (x),(v + h(6))6) + o(\\g + h(S)\\S) 
where lim^-^o = 0. By first dividing S and 


then taking the limit of S f 0 on both sides of the equa¬ 
tion above, we have 

(r 7 + h(5))5) - V(x)) = (V7Vj{x),g). (20) 

Since V satisfies Lipschitz condition, the limit in (20) co¬ 
incides with DV (x; g) and hence DV (x; g) = DVj (x; g). 
The above argument also holds for f, which completes 
the proof. □ 

Theorem 2 (PSCLF Theorem II) The closed-loop sys¬ 
tem (2) under the switching law v = <S[F] is globally 
asymptotically stable including sliding motions if for any 
x £ dfl n dV , there exists j £ J(x) such that 

DVj(x;f q (x)) < ~W(x),Vq £ Mj(x) nX“ m (x). (21) 

PROOF. Let x(t) be the closed-loop trajectory of sys¬ 
tem (2) under the switching law v = <S[V"]. We want 
to show that Dy(x(f);i(f)) < —W{x(t)) for t £ R + 
a.e.. By Theorem 1, we only need to show the inequality 
for the part of the trajectory that involves sliding mo¬ 
tions. Let (t i, £ 2 ) be the time interval that sliding motion 
occurs, i.e., x(t) = J2 ie x am { x {t)) oii(x(t))fi(x{t)) where 
a i{ x< f)) = 1 ,f £ (tiM) a-e.. The proof is 
divided into two cases: i) the ones on the smooth region 
({ Clj }j ) and ii) the ones on the nonsmooth boundary 
(5f2). We want to show that for both cases V satisfies 

DV(x(t)- : x(t)) < —W(x(tj) for t £ (< 1 X 2 ) a.e.. (22) 

To simplify notation, we will use x instead of x(t) in the 
rest of the proof. 

For case i), V is continuously differentiable at x 
by (11) and therefore Dk(x; g) = (VV(x),g) is 

affine with respect to g. As a result, DV(x;x) = 
Yji£i sm (: r ) a i( x )DV(x-, fi(x)). By the definition of l sm 
in (3) and the construction of v in (13), we have 
DV(x; fi(x)) = min ge g DV (x; f q (x)),Vi £ X sm (x). 
Then the desired condition (22) follows directly from 
the decreasing condition of V in (8). 

For case ii), According to condition (21), there exists 
some j* £ J(x) such that 

DVj. (x;/q(x)) < -W(x),Vq £ AT).(x) nX“ m (x). 

By part ii) of Lemma 4, 

DVj* (x; f q (x)) < —W(x),Vq £ Mj*(x). 

The above two inequalities give 

DVj* (x; /i(x)) < —W(x),Vi £l“ m (x). (23) 
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To stay on the nonsmooth boundary dfl, the velocity 
satisfies x £ Dje j(x) T x (dClj). By Lemma 6 , DV (x; x) = 
DVj * (x; x). As we know, DVj » (x; 77 ) is affine with respect 
to r), which gives 

DV(x;x)= ^ ai(x)DVj*(x;fi(x)) 

i(zX S m(x) 

= ^ ai(x)DVj* (x;/j(x)) <-W(x), 

where the second equality is due to the fact that a, (x) = 
0, Vi G T sm (x)\I“ m (x) and the last inequality follows 
from (23). Therefore, we complete the proof for case 
ii). □ 


According to Theorem 2, the switching law v = «S[V] 
can guarantee closed-loop stability including sliding mo¬ 
tions if V is a PSCLF and condition (21) is satisfied 
on dfl fl dV. For any given PSCLF, the switching law 
can be constructed accordingly and condition ( 21 ) can 
then be checked without requiring further runtime tra¬ 
jectory dependent information. Note that the index set 
A ij(x) PI X“ m (x) is a strict subset of Q and thus the 
condition ( 21 ) is less conservative than checking the in¬ 
equality for all the subsystems q £ Q. In fact, for many 
cases, the index set Mj(x) is empty for some j £ J(x), 
for which ( 21 ) holds trivially. 

We end this subsection by revisiting Example 3 to illus¬ 
trate how to check condition (21) when the switching 
surface dV partly coincides with the nonsmooth surface 
dLl. Recall that the nonsmooth surface of V is Si U S -2 
where Si = (x G R 2 : Xi - X 2 = 0} and S 2 = {x £ 
K 2 : Xi + 3x2 = 0}. Since I sm (x) = {1, 2}, Vx £ Si, Si 
is also a switching boundary, i.e., Si C dV. It can be 
easily verified that S 2 is not a switching boundary and 
therefore dfl fl dV = Si. By the relations summarized 
in (17), we have Ali(x) = {1}, Al 2 (x) = (2},Vx £ Si. 
If we fix j = 2, then the index set to be checked is 
Ad^x) nX“ m (x) = {1}. By the directional derivatives 
computed in (16), we have DV^x; /i(x)) = —6||x|| 2 and 
from which condition (21) is verified. Therefore, the sys¬ 
tem is globally asymptotically stable including sliding 
motions under the switching law v = S[V]. 

3-4 Important Special Cases 

This section discusses two important special cases of 
Theorem 2. The first one is when the CLF V is smooth 
over the entire state space, and the second one is when 
V is obtained by taking the pointwise minimum of a fi¬ 
nite number of smooth functions. Both cases have been 
studied in the literature. Using Theorem 2 we are able 
to obtain stronger results in a unified way. 


Corollary 1 (Smooth CLF) IfV is smooth and satisfies 
all the conditions (6), (7), and (8), then the closed-loop 
system (2) under the switching law v = «S[V] is globally 
asymptotically stable including sliding motions. 

When V is smooth, the nonsmooth boundary dfl is 
empty and condition (21) holds trivially. Thus, the 
above corollary follows immediately from Theorem 2. It 
is worth mentioning that with smooth V, the closed-loop 
vector field under v = <S[V] can still be discontinuous 
with trajectories involving sliding motions. Therefore, 
Corollary 1 is not a direct consequence of the classical 
CLF results. 

We next consider a special class of PSCLFs that are 
obtained by taking the pointwise minimum over a finite 
number of smooth functions. 

Definition 12 (Pointwise Minimum CLF) Let V pm j : 
R" — 5 - R +; j £ N m , be nonnegative continuously differ¬ 
entiable functions onR". The function defined by 

V P m(x) = min V pm j(x), x £ R”, (24) 

is called a pointwise minimum control-Lyapunov function 
(PMCLF) if 

(1) ftpmj 0, Mj £ N m where f } pm ,j — {x £ R n : 

Vprn.j (x) < ^pm,k{x),Wk f j \? 

(2) V prn satisfies conditions (6), (7), and (8). 

The first condition in the above definition ensures that 
every smooth function contributes nontrivially to 

the pointwise minimum. Due to the smoothness of V pm j , 
each set D pmj is open with continuously differentiable 
boundary <9Q pm . r Therefore, a PMCLF V pm is always 
piecewise smooth and it is also a PSCLF. 

Corollary 2 (PMCLF Theorem) IfV pm is a PMCLF, 
then the closed-loop system (2) under the switching law 
v = S[V pm \ is globally asymptotically stable including 
sliding motions. 


PROOF. Since tp m is piecewise smooth, it suffices to 
show that V pm satisfies condition (21). We pick an ar¬ 
bitrary point x £ dLl n dT>. Since I sm {x) 0, part 
ii) of Lemma 3 guarantees the existence of a j £ J (x) 
such that Mj(x ) nl sm (x) 7 ^ 0. We now fix this j and 
pick an arbitrary q £ Mj{x) PI l sm (x). If M){x) PI 
X sm (x) = 0, then condition (21) holds trivially. Now as¬ 
sume that ■Ad'j(x) Pll sm (x) 0 and pick an arbitrary 
q' £ Mj(x) PI I,™ (x). Then, it follows again from part 
ii) of Lemma 3 that there exists a j' £ J{x),j' j such 
that q' £ Mj'(x) Pi I S m (x). Next, we want to show that 
for x £ dfl, if i £ Mk{x) Pll“ m (x) for some k £ J(x), 
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then 

DV pni} k(x-, > DVp mij (x; G J(x). (25) 

By the regularity condition in Definition 5, there ex¬ 
ists e,; > 0 such that VS G (0, e^), x — Sfi(x ) G XV 
It follows from x G dfl k H dT>i that x — Sfi(x ) G 12^. 
By the definition of V pm in (24), V pm ^{x - $fi(x)) < 
V pm j(x — Sfi{x)),Vj G N m . According to the con¬ 
tinuity of Vpm* V p m,k(x) = Lpm, j(x), Vj G J {P). It 
follows that lim. 54,0 y(V pm ,fe(a;) - V pia , k (x - Sfi(x))) > 
lim^^o j(Vpm ,j(x) - V pm j(x - 5fi(x ))) and from which 
we proved (25). Apply (25) to some q' G Mj(x) nlf m (x) 
where q' G My{x), we have 

DVpuij' (x-, fq' (x)) > DVpmj(x; f q > (x)). (26) 

From part ii) of Lemma 4, DV pnl j'(x-, f q '(x)) < —W(x). 
Together with (26), DV pin j(x; f q >(x)) < —W(x), Vq' G 
Mj(x ) nXg m (x), which verifies condition (21) and there¬ 
fore completes the proof. □ 


Corollary 2 indicates that if a PMCLF is used, closed- 
loop stability including sliding motions can always be 
guaranteed without any extra condition on dT> 0 <912. 
Such a result holds for any switched nonlinear system 
and any PMCLF (not necessarily piecewise quadratic). 
It represents an important contribution on its own. 

4 Application Examples 

The PSCLF approach, along with the stability re¬ 
sults, provides a unified framework to design stabilizing 
switching laws with a systematic consideration of slid¬ 
ing motions. Once a PSCLF is found, the design of the 
switching law and the stability analysis of the closed- 
loop system follow directly from our PSCLF results. 
The search for a PSCLF can often be done numerically 
through proper parametrization of the PSCLF. Similar 
ideas have been studied extensively for switched linear 
systems (SLSs) with quadratic or piecewise quadratic 
CLFs [Ij, |2(j. In this section, we will first briefly show 
that the proposed framework can be used to recover and 
extend many existing methods for SLSs in a unified way, 
and then we will use a numerical example to illustrate 
its application in stabilization of switched nonlinear 
systems. 

4-1 Applications in Switched Linear Systems 

We first consider a general switched linear system (SLS) 
given by: 

x(t) = A a{t) x(t),cr(t) G Q = {!,•■• ,M},teR+, (27) 


where a : R + —>• Q denotes the switching signal and 
{Ahe q are constant matrices. Note that asymptotic 
stability is equivalent to exponential stability for SLSs 
[l8j |. This fact is implicitly used in some parts of the 
following discussions. 

V 1.1 Quadratic Switching Stabilization 

A well studied stabilization problem for SLSs is the so- 
called quadratic stabilization problem. System (27) is 
called quadratically stabilizable if there exists a switching 
law v under which the closed-loop system has a quadratic 
Lyapunov function (2lll25l |. Using our framework, we call 
a SLS quadratically stabilizable if it admits a quadratic 
CLF of the form V{x) = x T Px , x G K ra . In this special 
case, it can be easily verified that the PSCLF conditions 
in (6), (7), and (8) are equivalent to 

Py 0, and mmx T (Af P + PAi)x < 0,Vx 0. (28) 
ieQ 

This coincides with the strict completeness condition 
proposed in (25}. It can be easily verified that the follow¬ 
ing condition is a sufficient condition to ensure (28): 

P y 0, and ^ caAi'j P + p(^^ cqA^ -< 0, 
ieQ ieQ 

for some a,/ > 0 with a, = 1. (29) 

ieQ 

The above condition is known as the “stable convex 
combination” condition 0, [34} and can be used to find 
P by solving a linear matrix inequality (LMI) feasibil¬ 
ity problem. Once a quadratic function satisfying (28) 
is found, the stabilizing switching law can be obtained 
immediately using our framework. Since the function 
V(x) = x T Px is globally smooth, the switching law de¬ 
sign reduces to v{x) = argmin^Q x T (Af P + PAi)x. 
Corollary 1 guarantees the closed-loop system under 
this v is asymptotically stable including sliding motions. 
Therefore, our framework can be used to recover most of 
the existing results in quadratic switching stabilization 
problems. 

4-1-2 Piecewise Quadratic Switching Stabilization 

As a natural extension of quadratic stabilization, we 
can consider piecewise quadratic functions as can¬ 
didate CLFs. A well known result along this direc¬ 
tion is the largest-region switching strategy [U, 122} . 
whose construction depends on two key components. 
The first one is a collection of regions defined by 
X2j = {iG R" : x T Hi,x >0},tG Q. The second compo¬ 
nent is a collection of quadratic Lyapunov-like functions 
Vi{x) = x T Pi.x, x G 1", i G Q. Note that for each i G <2, 
the matrices Hi and Pi are symmetric but may not be 
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positive or negative semidefinite. Given the two compo¬ 
nents, the largest region switching strategy is defined by 
v(x ) = argmax^Q x T HiX. It has been shown that this 
switching law guarantees closed-loop stability excluding 
sliding motions under the following four conditions. 

(Ql’) The union of the regions covers the entire space, 
i.e., Uiggfij = M”; 

(Q2’) Vi is positive definite on Clj; 

(Q3’) The Lie derivative of Vi along the vector field of 
subsystem i is negative definite on fb; 

(Q4’) Vi(x) = Vj{x) on {a; £ : x T HiX = x T Hjx}, for 

all i,j £ Q. 

Using the 5-procedure Q, the matrices { H,. Pi} iG Q and 
hence the largest-region switching strategy can be found 
by solving some bilinear matrix inequalities (BMIs). The 
derivation of these BMIs can be found in [2l], [22|. 

The largest-region switching strategy mentioned 
above can be quite conservative. First, the number of 
Lyapunov-like functions has to be equal to the number of 
subsystems, which is a restrictive assumption. Second, 
the decreasing condition for V, in region Clj is also conser¬ 
vative. Note that the regions {Sli}j e g are not mutually 
exclusive and they differ from the actual switching re¬ 
gions Clj = {a; £ K” : x T HiX > x T HjX,Vj / i},i 6 Q. 
Therefore, requiring V to decrease on f \ would be a 
better choice. Last, selecting subsystem based on the re¬ 
gion matrices, although can simplify stability analysis, 
is less effective. Roughly speaking, the switching control 
should be chosen to decrease Lyapunov-like functions 
to achieve a better closed-loop stability performance. 

Our framework can be used to tackle these issues and ex¬ 
tend the largest-region method in a systematic way. For 
example, we can consider a piecewise quadratic function 
V : R" —>• R + with m partitions defined by Clj = {x £ 
R n : x T HjX > x T HkX,\/k ^ j,k £ N m },j £ N m . The 
restriction of V to Clj is assumed to take a quadratic 
form Vj (x) = x T PjX, x G Clj. It can be verified that 
this function V will be a PSCLF for system (27) if the 
following conditions hold. 

(Ql) For each j G N m , Vj is positive definite on Clj-, 

(Q2) For each j G N m , min ie g x T (AfPj + PjAi)x < 
0,Vx G Clj; 

(Q3) Vk(x) = Vj(x), Vx G dClk fl dClj, for all k,j G N m . 

When V is a PSCLF, Theorem 1 guarantees that the 
closed-loop system under the switching law v = 5[U] is 
asymptotically stable excluding sliding motions. Stable 
sliding motion can be also guaranteed by an additional 
condition 

(Q4) Vj G N m , 3k G N m such that x T (Af Pk + PkAi)x < 
0, Vx G dClj, Vi G Q. 


The above condition (Q4) guarantees (21) for the 
PSCLF V as required by Theorem 2. A set of BMIs can 
be derived using the 5-procedure to guarantee these 
four conditions. 

Theorem 3 If there exists real, symmetric matrices 
Pj,Hj,j G N m , and real numbers£jk, Kjkt G R, r]i, r] 2 > 
0? 'Yjk, Pjk, Cijk Cl 0 7 Q-ij G [0,1], j, k, t G G Q such 

that Y^ieQ a ij = 1 yj G N m , and 

Pj~i111 h T, 7 jk{Hj - Hk),Vj € N to ; (30a) 

fce N m 

( Pj + Pj ( + V2Pj 

ieQ ieQ ( 30b ) 

^ Pjk(Hk - Hj),\/j G N m ; 
keNm 

Pj = Pk + £jk(Hj - H k ),Vj, k G N m ; (30c) 

Qjk(Aj Pk + PkAi)+ 

(30d) 

A ijkt{Hj — H t ) -< 0,Vz G Q, Vj G N m , 

k,te N m 

then under the switching law v = 5[U], system (27) is 
globally asymptotically stable including sliding motions. 


PROOF. By the 5-procedure, it is easy to verify that 
(30a) and (30c) imply condition (Ql) and (Q3), respec¬ 
tively. For any x G Clj, we have x T HjX > x T HkX,Mk G 
N m and hence j3jkX T {Hk — Hj)x < 0 for any (3jk > 0. 
By left multiplying x T and right multiplying x on both 
sides of (30b), we have 

DVjfx;^ aijAiX ^ +r] 2 Vj(x) < 0,Va; G Clj. 

iGQ 

Since J2ieQ = 1, it follows that for each j G N m , 

min DVj(x; Aix) < DVj (x; QyAjij 

ieQ 

< —ri2Vj(x),Vx G Clj, 

which verifies condition (Q2). Next, we want to show 
that (30d) implies condition (Q4). Since Qjk > 0, 

Cijk(Af Pk + PkAi — bijkt{Hj — H t )) -< 0 

ken™. teN m 

implies that there exists at least one k G N m such that 
AfP k + P k Ai - b ijkt(Hj - H t ) -< 0, otherwise 

the above inequality changes direction. We now fix this 
k and by the 5-procedure we know that x T (AfPk + 
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PkAi)x < 0,Vx £ d£lj, which verifies condition (Q4). 
Let A ijkt — ~Cijkbijkt j we obtain the form in (30d). □ 


The search for a piecewise quadratic CLF can thus be 
formulated as a feasibility problem for the BMIs given 
in Theorem 3. Such a result is more general than the 
largest-region switching approach [22[ as it allows the 
number of the switching regions to be different from the 
number of subsystems, and guarantees closed-loop sta¬ 
bility including sliding motions. 


4-1.3 Switching Stabilization with Composite Control- 
Lyapunov Functions 


'IjkiPjkiCijk — 0? otij £ [0, l],j, £ Q such 

that A A v- q ( x-ij — 1, Vj £ N m , and 


Pj ~ T)il >Z ^2 7jk{Pj - Pk),Vj e N m ; 

fceN m 

( P i + Pj ( y~l + mPj 

ieQ ieC 

A G N m ; 

feGNm 

c«fc(4’ft+^)+ 

fceNm 

A ijkt(Pj — Pt) A 0, Vi £ <2, Vj £ N m 

fc,ieNm 


(31a) 


(31b) 


(31c) 


Another class of CLFs that have been studied for SLSs is 
composite quadratic functions, which are defined by tak¬ 
ing the pointwise minimum, pointwise maximum, and 
convex hull of a finite number of quadratic functions [3l ■ 
The stabilization results based on these three classes of 
composite CLFs are also special cases of our PSCLF 
framework. Moreover, they can be further extended and 
strengthened using our framework. As an example, we 
consider a pointwise maximum CLF defined by: 

14m — max Vj(x), x £ R", 
jSN m J 

where Vj(x) = x T PjX. Obviously, 14m is a piecewise 
quadratic function with to partitions defined by f lj = 
{x £ R n : x T PjX > x T PkX,Mk ^ j,k £ N m }, j £ N m . 
It can be verified that 14m will be a PSCLF for system 
(27) if the following conditions hold. 

(Ml) For each j £ N m , Vj is positive definite on fl,-; 

(M2) For each j £ N m , min ie QX T (AfPj + PjAt)x < 
0, Vr £ ttj. 

When P p m is a PSCLF, Theorem 1 guarantees that the 
closed-loop system under the switching law v = S [14m] 
is asymptotically stable excluding sliding motions. Sta¬ 
ble sliding motion can be guaranteed by an additional 
condition 

(M3) Vj £ N m , 3k £ N m such that x T (Af P% + PkAi)x < 
0, Va; £ dttj,\/i £ Q. 

Note that the above condition implies condition (21). 
It requires x £ dflj rather than x £ dflj fl dT>i for 
ease of BMI derivation. The construction of P p m and the 
corresponding switching law design can be formulated 
as a BMI feasibility problem as stated in the following 
theorem. 


then under the switching law v = S [14m]. system (27) is 
globally asymptotically stable including sliding motions. 

PROOF. By the 5-procedure, condition (Ml) is guar¬ 
anteed by (31a). Similar as the argument of (30b) =>■ 
condition (Q2) shown in the proof of Theorem 3, we 
have (31b) implies condition (M2). We are left to show 
that (31c) implies condition (M3). Since (jjfc > 0, 
SfceN m (ijk(Af Pk + PkAi — bijkt(Pj Pt)) A 0 

implies that there exists at least one k £ N m such that 
Af Pk + PkAi — ^2 bijkt(Pj ~ Pt) 0. (32) 

teN m 

We now fix this k. Apply the 5-procedure to (32), we 
have x T {AjPk + PkA^x < 0 on {x £ R n : x^PjX = 
x T P t x,t £ N m } = U tS N m (i9f2j D dClt)- Due to the fact 
that dQj = UtgN m (c>f2j n 9fl t ), condition (M3) is ver¬ 
ified. Let A ijkt — —Cijkbijkt , we obtain the form in 
(31c). □ 

The authors in pjj also derived a set of conditions to 
guarantee closed-loop stability including sliding mo¬ 
tions. The conditions are mostly the same as ours except 
for the one on the nonsmooth boundary dfl. The BMI 
conditions in 0 are derived only for the case where 
P p m is composed from two quadratic functions (i.e. 
to = 2), which is more restrictive than condition (31c) 
obtained using our framework. It is worth mentioning 
that similar extensions can be made for the case with 
convex-hull composite CLFs. As for the pointwise min¬ 
imum CLF case, the result in []j[ is a very special case 
of our general result given in Corollary 2. All these dis¬ 
cussions further illustrate the importance and unified 
nature of the proposed PSCLF framework. 

4-2 A Numerical Nonlinear Example 


Theorem 4 If there exists real, symmetric matrices Switching stabilization of switched nonlinear systems 
Pj,j £ N m , and real numbers Xijkt £ ? 7 i, 772 > 0, has not been adequately studied in the literature. Here, 
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(a) (b) (c) 




(d) (e) 

Fig. 4. (a): The phase portrait of subsystem 1, (b): The phase portrait of subsystem 2, (c): The phase portrait of the closed-loop 
system under the switching law v, (d): The closed-loop trajectory of system (33) with xo = [1,0] T , (e): The closed-loop 
trajectory of system (33) with xo = [0, — 1] T . 


we use a numerical example to illustrate the applica¬ 
tion of our framework in this area. Consider the fol¬ 
lowing switched nonlinear system with two subsystems, 
x[t) = fi(x(t)),i = 1,2, where 


h(x) 


—*i — x 2 
*f + 0.5x2 


J2(x) 


*1 — x 2 
x\ — 1.5*2 


(33) 


For this simple switched polynomial system with sub¬ 
system vector fields shown in Fig. 4(a) and Fig. 4(b), we 
consider a smooth polynomial CLF:!/(*) = xf + 2x\ 
under which the closed-loop vector field is shown in 
Fig. 4(c). Apparently, V satisfies the positive definite 
condition and the radially unbounded condition given 
in (6) and (7). In addition, we have DV(x;fi(x)) = 
— (4*i - 4*|) - 2*|, DV(x; f 2 (x)) = (4*f - 4*|) - 2x\. 
The above implies min^ DV{x\ fi{x)) < — 2*|, which 
verifies the decreasing condition in (8). Therefore, V 
is a special case of the proposed PSCLFs. Since it is 
smooth, the corresponding switching law in (14) takes 


the following form: v(x) = argmin^ DV(x; fi(x)). By 
Corollary 1, we can conclude that the closed-loop sys¬ 
tem * = / y ( x ) (*) is asymptotically stable including slid¬ 
ing motions. The result is also illustrated through sim¬ 
ulations. In particular, Fig. 4(d) and Fig. 4(e) show the 
closed-loop trajectories starting from * 0 = [1,0] T and 
*o = [0, — 1] T , respectively. Both trajectories exhibit 
sliding motion behavior, which can be observed in the 
zoom-in box. We introduce hysteresis band in the sim¬ 
ulation to deal with the discontinuous closed-loop vec¬ 
tor field. Both trajectories converge to the origin under 
sliding motion as expected. 

5 Conclusion 

In this paper, we proposed a piecewise smooth control- 
Lyapunov function (PSCLF) framework to study switch¬ 
ing stabilization problems. We formally introduced the 
concept of stability including or excluding sliding mo¬ 
tions and sufficient conditions in terms of PSCLFs were 
derived for the two stability notions. A constructive way 
to design a stabilizing switching law based on a given 
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PSCLF was also developed. We showed that such a con¬ 
trol law can guarantee the closed-loop stability exclud¬ 
ing sliding motions, and it can also ensure closed-loop 
stability including sliding motions under an additional 
condition on the nonsmooth surface of the PSCLF. We 
also showed that for smooth CLFs and pointwise mini¬ 
mum CLFs, stable sliding motions are always guaranteed 
without any additional condition. The proposed frame¬ 
work can be used to obtain and extend most of the ex¬ 
isting results in stabilization of switched linear systems. 
In addition, it provides a systematic way to study stabi¬ 
lization of switched nonlinear systems. Future research 
will focus on necessary stabilizability conditions and con¬ 
verse control-Lyapunov function theorems for switched 
linear systems. 
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